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■ Abstract 

ft Diquarks with J p = 0,1 containing a heavy (charm or bottom) quark and a light quark are investigated 

using QCD Laplace sum rules. Masses are determined using appropriately constructed gauge invariant correlation 
functions, including for the first time next-to-leading order perturbative contributions. The J F = + and 1 + 
charm-light diquark masses are respectively found to be 1.86 ± 0.05 GeV and 1.87 ± 0.10 GeV, while those of the 
+ and 1 + bottom-light diquarks are both determined to be 5.08±0.04GeV. The sum rules derived for heavy-light 
diquarks with negative parity are poorly behaved and do not permit unambiguous mass predictions, in agreement 

i-ft with previous results for negative parity light diquarks. The scalar and axial vector heavy-light diquark masses 

are degenerate within uncertainty, as expected by heavy quark symmetry considerations. Furthermore, these 
mass predictions are in good agreement with masses extracted in constituent diquark models of the tetraquark 
candidates X(3872) and lb(10890). Thus these results provide QCD support for the interpretation of the X(3872) 

_^ and Vb(10890) as J = 1 ++ tetraquark states composed of diquark clusters. Further implications for tetraquarks 

among the heavy quarkonium-like XYZ states are discussed. 

v]q 1 Introduction 

The discovery of the X(3872) by the Belle collaboration ft] and its subsequent confirmation by the CDF [SJ, D0 [3J, 
BABAR 4 and LHCb [5] collaborations initiated a new era in hadron spectroscopy. Since then, hadrons have been 
^+ found in the charmonium and bottomonium spectra that are difficult to reconcile as conventional heavy quarkonia. 

These are called heavy quarkonium-like or XYZ states, and a comprehensive review of the current experimental 
situation is given in Ref. [6 . The X(3872) exemplifies the difficulties in interpreting these states: its mass is 
M = 3871.68 ± 0.17 MeV, its width is T < 1.2 MeV 6^ and the LHCb collaboration has clearly established that its 
^ quantum numbers are J PC = 1 ++ 7J. These properties pose problems for a conventional charmonium interpretation 
of the X(3872) [8 . Given the proximity of its mass to that of DD* , the X(3872) has been widely interpreted as a 
four-quark molecular state [5J [TUJ QTJ Q21 1131 [HI HS1 H] • A complementary interpretation is that the X(3872) is a 
tetraquark [T71 [TS1 [THl HOI EI] ■ In addition to the X(3872), several XYZ states that are four-quark candidates are 
discussed in Ref. [2"2"| . 

Molecules and tetraquarks have very different internal quark structures. In the molecular scenario, two color- 
singlet mesons form a weakly bound conglomerate, whereas in the tetraquark scenario a diquark and anti-diquark 
form a tightly bound four-quark state. A diquark is a strongly correlated pair of quarks within a hadron (see Ref. [23] 
for a review of applications). Because single gluon exchange leads to an attractive interaction between quarks in 
a color anti-triplet configuration, diquarks are identical to anti-quarks in terms of color. In Ref. |24j all possible 
diquark configurations were classified and it was shown that due to spin interactions, the scalar is the most strongly 
bound, followed by the vector. However, these spin interactions scale as the inverse of the quark mass, and hence 
scalar and vector diquarks that contain one or more heavy quarks should be degenerate. 

The tetraquark and molecular currents used in QCD sum rule analyses are related through Fierz transformations, 
leading to ambiguities in their interpretation which can be addressed through the diquark scenario [25) (Ref. [26] 
provides a review of the numerous QCD sum rule studies of tetraquarks and molecules among the XYZ states). 
In addition, the renormalization of four-quark operators is complicated by operator mixing [57, 28 . Conversely, 
the renormalization of the diquark operator is multiplicative and has been studied to two-loop order [29 . For this 
reason, QCD sum rule studies of diquarks can be extended to higher orders much more easily. The first QCD sum 
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rule studies of diquarks were given in Refs. [3UJ ESJ, followed by Refs. [351 1321 f33J. The Bethe-Salpeter [371 135) . 
Dyson-Schwinger [39) . and effective field theory approaches [30] have also been used to determine diquark masses. 
Ref. [33] used QCD sum rules to investigate heavy- light diquarks with J p = + and f + . In this paper we will build 
upon previous work by including next-to-leading order perturbative contributions and negative parity diquarks in 
our analysis. 

Diquarks are clearly not hadrons, thus their masses must be regarded as constituent masses. Constituent diquark 
models have been used to study tetraquarks among the XYZ states. In Ref. 17 Maiani et al. interpret the X(3872) 
as a tetraquark composed of charm-light diquarks, and using its mass determine both the scalar and vector charm- 
light constituent diquark masses to be 1.93 GeV. Ref. [41] points out that the Zf(3895), which was very recently 
discovered by the BESIII [32] collaboration and quickly confirmed by the Belle [33] and CLEO [33] collaborations, 
was predicted in Ref. |17j . The confirmation of this charged charmonium-like state strongly supports the existence 
of hadrons outside the constituent quark model. Similarly, Ali et al. [33] interpret the Y"b(10890) discovered by 
Belle [35] as a tetraquark composed of bottom-light diquarks, determining the scalar and vector bottom-light diquark 
masses to be 5.20 GeV. The charged bottomonium-like states Z^ (10610) and Z^ (10650) [36] are also suggested to be 
tetraquarks. Important features of the analyses in Refs. [T7J133] are the use of heavy-light diquarks whose constituent 
masses are extracted from fits to tetraquark candidates and the equality of scalar and vector heavy-light diquark 
masses. In this paper we seek to determine if these heavy-light diquark masses are supported by QCD sum rule 
analyses, thereby providing a QCD-based test of the heavy-light diquark model of tetraquark states. Because our aim 
is to compare our results with the heavy- light diquark masses determined in Refs. [T7j[34], our focus is on heavy-light 
diquarks. 

The remainder of the paper is organized as follows: in Section [2] we calculate the J p = 0* , l* heavy-light 
diquark correlation functions, in Section [3] we construct and analyze the corresponding QCD Laplace sum rules, and 
in Section [4j we make concluding remarks and discuss the phenomenological implications of our results. 

2 Heavy-Light Diquark Correlation Function 

The heavy-light diquark correlation function is defined as 

n(Q 2 ) = ifd A xe i ^{Q\T[J a (x)S auJ [x,Q] 4 (0) ] |0) , (1) 

where Q 2 = —q 2 is the Euclidean momentum, and a, u> are color indices. The heavy-light diquark currents are 

J a — ^afiiQp,COq 1 , (2) 

where C is the charge conjugation operator, T denotes the transpose, Q is a heavy (charm or bottom) quark field, and 
q is a light quark field J3UJ EJ. The Lorentz structures O = 75 , 1 , 7^ , 7^75 respectively couple to scalar ( J p — + ) , 
pseudoscalar (0~), axial vector (1 + ), and vector (1~) heavy-light diquarks. We denote these as S, P, A, and V, 
respectively. The axial vector and vector correlation functions are given by 

II (A < V) (q) = ^ L l { t f- r) lC V) (q) , (3) 

where d is the number of spacetime dimensions. Following Refs. [301 [3TJ [25] [33l [32] ? the diquark correlation function (1IJ) 
includes a path-ordered exponential, or Schwinger string, defined as 



S aiLl [x , 0] = P exp 



ig^f- J* dz» Al(z) 



(4) 



where P denotes path-ordering and g = y/Aira is the strong coupling. Ref. [311] demonstrated that the correlation 
function (fTl) is gauge invariant to leading order for light quark currents. We will show that this is also true for 
heavy-light diquark currents (pi. 

First we calculate perturbative contributions to the heavy-light diquark correlation function, which are shown in 
Fig. [I] We include 0(a) perturbative contributions that have not been calculated previously. To leading order the 
Schwinger string Q generates a trace over the color indices in (ITJ , and this trace has been performed in calculating 
perturbative contributions. We will also consider a higher order contribution from the Schwinger string that is 
gauge dependent and should cancel the gauge dependence of the perturbative contributions. Thus perturbative 



Figure 1: Feynman diagrams representing the leading order and next-to-leading order perturbative 
contributions to the heavy- light diquark correlation function (TlJ> . An insertion of the diquark 
current is represented by the <8> symbol, bold lines represent heavy quark propagators, thin lines 
represent light quark propagators, and wavy lines represent gluon propagators. An additional 
diagram where the light and heavy quark propagators are exchanged is not shown. These and all 
subsequent Feynman diagrams were created using JaxoDraw [45] . 



contributions are calculated in a general covariant gauge in order to verify the gauge independence of the correlation 
function ([I]). The gluon propagator is taken to be 
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where a denotes the gauge parameter and the functions D^J , D^J are defined for later convenience. As in Refs. [461147] 
we calculate the entire correlation function, rather than only the imaginary part. This approach is essential in order to 
deal with gauge invariance and renormalization issues properly in this calculation. Results for the loop integrals that 
are encountered are given in Refs. [1H1 SHI EI] • The number of distinct integrals to be calculated can be significantly 
reduced using the Mathematica package Tarcer |53j , which implements the generalized recurrence relations developed 
in Refs. 54, 55 . Finally, the epsilon expansion can be performed using the Mathematica package HypExp [5H1I57J . 
Using the MS scheme and working in d — 4 + 2e dimensions, the perturbative result for each channel can be 
parametrized as 
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Here the subscript B indicates bare quantities, i = S , P , A , V denotes each distinct channel, Li 3 and Li 2 denote 
the trilogarithm and dilogarithm functions 58], and we have omitted terms corresponding to dispersion relation 
subtraction constants. The coefficients bi are functions of to which are given for each channel in Table [T] 

Some comments must be made regarding the form of (pf. First, terms proportional to the gauge parameter a 
have been retained to allow comparison with contributions from the path-ordered exponential Q, so as to ensure 
that the correlation function (fl]) is gauge invariant. Second, the term 64 in ([6]) is a non-local divergence that cannot 
be removed through application of the Borel transform when the sum rules are constructed. This term must be dealt 
with through renormalization, necessitating inclusion of the terms 61, 6 2 , and 63 which will lead to renormalization- 
induced contributions. A similar methodology was also needed in Ref. [S5]. Finally, note that (pi) has a branch cut on 
to € (—00 , — 1], as it must. However, after using the package HypExp some functions are generated that do not have 
this branch structure. This anomalous branch structure is spurious and is eliminated when polylogarithm identities 
are used [58] , 

We now turn our attention to contributions from the Schwinger string Q . Following Ref. [30] [3Tj , we define 
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Table 1: Coefficient functions bt for the bare perturbative result (KM). Here L m = log ^ • 

where : : denotes normal ordering. As in [301131] the integration path between points and x in Q has been chosen 
to be a straight lineF] As mentioned earlier, the leading order term in ^f\ leads to a trace over the diquark current 
color indices in (fTl), which was done in calculating (km). To the order that we are working, the quadratic term in g 
is irrelevant because it cannot be used to form a gluon propagator. However, the linear term leads to a non-trivial 
contribution to the correlation function, which is shown in Fig. [2j This contribution has the form 



^sL s (Q 2 ) ~ J d 4 x e*** J d A z J 1 d£ x* Dtfifr - z) 



(8) 



where q is the external momentum, z denotes the location of the quark-gluon interaction in Fig. [21 and the ellipses 
indicate ^-independent terms that are not shown. The £ integration in Q8J) cannot be evaluated readily. In momentum 
space the gluon propagator in ([8]) unavoidably leads to terms of the form 

1 In Rcf. 31 it was argued that any deviations from a straight line would correspond to additional Wilson loops, and hence would not 
correspond to the lowest energy configuration. 




Figure 2: Feynman diagrams representing the contribution of the Schwinger string to the heavy- 
light diquark correlation function. An insertion of the Schwinger string operator is represented by 
the © symbol. The dashed line is not a particle propagator. Instead, it indicates the straight line 
integration path between points and x used in equation (TF|. An additional diagram where the 
light and heavy quark propagators are exchanged is not shown. All other notations are identical 
to Fig. [1] 

where the £ and loop integrations are coupled. Most of these integrals can be decoupled using the scaling properties 
of (i-dimcnsional momentum integrals [50] . but unfortunately a few cannot be. However, for the gauge dependent 
terms in Ety this obstacle can be circumvented. Note that the quark propagator in (p)]) satisfies the identity 

h ■-^-S(q-k 2 - fa) = -^S(q-k 2 - £ki) . (10) 

The DftJ^ki) part of the gluon propagator (l5| provides a factor of fc^, and hence the £ integration in fcty can be 
performed using (10). Note that this approach cannot be used to calculate terms in ^ that correspond to the 



DfiJiki) piece of the gluon propagator ^. Based upon the result of Ref. [3DJ we have assumed that there is no 
contribution from the Schwinger string in Landau gauge. The remaining loop integrations can be performed using 
the same methods that were used to calculate the perturbative contributions. It should be noted that because of 
the bosonic nature of the diquark currents and gauge field, the integration over £ in fcfh must be symmetric about 
the point £ = oi meaning that the gauge configurations corresponding to £ and 1 — £ are equivalent. Thus the £ 
integration double counts and we have introduced an overall factor of | accordingly. 

In order to check the validity of the methods described, we have used them to reproduce the result of Ref. [3DJ , 
verifying that the correlation function ([lj for light diquark currents is gauge independent to order a. Using the 
approach described above we have calculated the gauge dependent contributions of the Schwinger string ufh to the 
heavy- light diquark correlation function (nl). We find that these precisely cancel the gauge dependent terms 612, 
&13 > &14 > an d &15 in the perturbative contribution (p)]). This verification of gauge independence emerges from the 
manifestly gauge invariant formalism of the Schwinger string, confirming that the heavy-light diquark correlation 
function (IT]) is suitable for use in a QCD sum rule analysis. 

Now we must renormalize the bare result (km). To the order that we are working, this can be done through 
renormalization of the heavy quark mass and the diquark current. The one- loop expression for the renormalized 
quark mass is [5T 

a 

m B = Z m m, Z m = H . (11) 

ire 

The renormalization of the scalar diquark current was studied in Ref. 29 . A distinct benefit of using a diquark current 
rather than a four quark current is that, unlike four quark currents, the diquark current renormalizes multiplicatively. 
In Ref. |29j it was shown that the renormalization factors of the scalar diquark and meson operators are proportional 
at one-loop level. This relationship can be extended to the pseudoscalar, axial vector, and vector channels in 
order to determine the renormalization factors of those diquark operators. Given our explicit demonstration of 
gauge independence and that the Schwinger string contributions are zero in Landau gauge [3D], we calculate the 
renormalization factors in Landau gaugejj The results are as follows: 

\m*=z$>\m B ', z P=^+^ z i p) = l+ i' Z « (A)=1 < z -°°= 1 - (12) 

Note that the axial vector and vector diquark operator renormalization factors are trivial, in analogy with the 
corresponding meson operators. Finally, the renormalized perturbative result for each distinct heavy-diquark channel 

2 This is the approach that was implicitly used in Rets. [30, 31 . 
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The heavy quark mass and strong coupling are implicitly functions of the renormalization scale fi, and the coefficients 
Cj are functions of w that are given in Table [2J 
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Table 2: Coefficient functions c* for the renormalized perturbative result ( |13[ ). All notations are 
identical to those in Table IT] 



The imaginary part of ( 13 ) can be easily determined via analytic continuation. The result is 
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where the coefficients fi are functions of x as given in Table [3] 
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Figure 3: Feynman diagrams representing the dimension-four quark condensate m q {qq) contribu- 
tion to the heavy-light diquark correlation function. Solid dots represent field condensates. All 
other notations are identical to Fig. [T| 

Now we consider contributions to the heavy-light diquark correlation function from the QCD condensates. Follow- 
ing Ref. [30, 51] . we calculate these contributions using fixed-point gauge techniques because the Schwinger string frf\ 
does not contribute to the condensates due to the x'M^ = gauge condition. We note that the manifestly gauge 
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Table 3: Coefficient functions fc for the imaginary part of the renormalized perturbative re- 
sult Ety. 



invariant nature of the correlation function (TTj) containing the Schwinger string implies that the fixed-point gauge 
results will be equivalent to those obtained in other methods |52) . First we consider the contribution from the quark 
condensate (qq) , which is shown in Fig. [3] For this contribution we find 
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(15) 



When the QCD Laplace sum rules are constructed in Sectional we will need to calculate the Borel transform B of ( 15 ) 
and all additional condensate contributions. The following result is useful in order to calculate Borel transforms of 
the condensate contributions [ST] 
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This result can be extended to cases where the denominator is raised to a higher power by differentiating ( 16 ) with 
respect to to 2 . Using this result, the quark condensate contributions to the sum rules are given by 
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Figure 4: Feynman diagrams representing the dimension-four gluon condensate (aG 2 ) contribution 
to the heavy-light diquark correlation function. An additional diagram where the light and heavy 
quark propagators are exchanged is not shown. All notations are identical to those in Fig. |3| 

Next, we determine contributions from the gluon condensate (aG 2 ) = (aGp„G£ v ) , which are shown in Fig. W| 
For these contributions we find 
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The Borel transforms of these are 
*G 2 
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In calculating (191 for the axial vector and vector channels we have not included the logarithmic term in (18). This 



term will lead to an imaginary part and hence the gluon condensate will have a continuum contribution in these 



channels. This can be calculated an identical fashion to (|14|), with the result 
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Figure 5: Feynman diagram representing one of the dimension-five mixed condensate (gqaGq) 
contributions to heavy-light diquark correlation function. All notations are identical to Fig. [T| 



The contributions of the mixed condensate (qcrGq) = (gq^-a^G^q) are 
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Note that (21 1 includes a term that arises from the fixed-point gauge expansion of the vacuum expectation value 
(q(x)q(0)) in Fig. [3] This is separate and distinct from the term that is represented in Fig. [5J The contributions of 
the mixed condensate to the sum rules can be calculated using (|16[), yielding 
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Finally we consider the dimension-six quark condensate, a(qq} 2 7 which arises purely from a higher order term in 
the fixed-point gauge expansion of the vacuum expectation value (q(x)q(O)) in Fig. pj For this we find 



n w - (o 2 

±x qqqq V^ 



16vr 

In 



a{qq)' 



(Q 2 + m 2 y 



u {F) - (o 2 ) = u {s) - (o 2 ) 

11 qqqq \^ I lx qqqq \^ J 



(23) 



,2\ _ _ TT (S) 



n^ (o 2 ) = -m bj - (o 2 ) if '- (o 2 ) = -u^'- (o 

qqqq V^ / qqqq v^ / ' qqqq V^ / qqqq v^ 



(V) 



t(S) 



where we have assumed vacuum saturation. The contributions of the dimension-six quark condensate to the sum 
rules are given by 
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We do not consider the dimension-six gluon condensate in this analysis. In Section [3] we will see that the gluon 
condensate is a sub-leading contribution to the heavy-light diquark sum rules, hence we expect higher-dimensional 
gluon condensates are suppressed and can be ignored. 



3 QCD Laplace Sum-Rule Analysis 

We now proceed to the QCD Laplace sum rules analysis of J p = ± , 1 heavy-light diquarks. Refs. [501 EI] are the 
original papers presenting the QCD sum rules technique, and reviews of its methodology are given in Refs. [62(163]. 
Using a resonance plus continuum model for the hadronic spectral function 



p^(t)=p lcs (t) + 9(t- SQ )ImIL(t), 
where sq is the continuum threshold, the Laplace sum rules take the form 
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where to is the hadronic threshold. The left hand side of (26 1 is given by 
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We now construct the heavy-light diquark sum rules. Using the results obtained above for the perturbative ( 14 1 



quark condensate (17), gluon condensate (19) and (20), mixed condensate (22), and dimension-six quark conden 



sate (24) contributions, the QCD Laplace sum rules are given by 

to 2 / >So/ "" 



*f(^o) 



7T 



(to z) 



Imll 



(i) 
port 



imn 



(i) 

GG 



z dz 



(28) 



■ Bg (k , r) + B« G (k,r)+ B% (k,r) + B^ q (k , r) 



The mass and coupling in ([28]) are implicitly functions of the renormalization scale p in the MS-scheme and renor- 
malization group improvement may be implemented by setting p = \j\fr [64) . In order to extract mass predictions 
for heavy-light diquarks we utilize a single narrow resonance model 



1 



p ICS (t) = f 2 6(t-M 2 ) 



(29) 



Eqn. ( 26 ) then yields 

n k (T,s Q ) = f 2 M 2k exv(-M 2 r), 
from which the heavy-light diquark mass M can be determined via the ratio 



M 



T^o (t, s ) 



(30) 



(31) 



Prior to extracting mass predictions we must discuss the QCD parameters occurring in the sum rules. We use 
one-loop MS expressions for the running coupling, charm and bottom quark masses: 



ot(p) 



:(M) 



l + A 



*(M) 



log 



uo 



, TO,(/i) = to ( — —— ) , m = m (p = to) 
\a(m)J 



In the charm-light diquark analysis we take 

M = M r = 1.77GeV, a(M T ) = 0.33 ± 0.01 , A = A c 
while in the bottom-light diquark analysis we use 



25 
12 



1.28 ± 0.03 GeV, 



23 



M = M Z = 91.188 GeV, a(M z ) = 0.1184 ± 0.0007, A = A b = — , fn b = 4.18 ± 0.03 GeV. 



(32) 



(33) 



(34) 



All of these parameters are taken from Rcf. [BJ, apart from A c and Aj, which are given in Ref. [51] . We set p = l/y/r 
in order to implement renormalization group improvement as described above. 
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We now specify the values used for the QCD condensates. Beginning with the quark condensate, we define 
to (2 GeV) 



m(qq) 



i q (2 GeV) 



m q {qq) 



(35) 



where m denotes the charm or bottom quark mass and we use the PCAC relation m q (qq) — — ^f 2 m 2 . The numerical 
values are again taken from Ref. [BJ: 



m q (2 GeV) = - [m u (2 GeV) + m d (2 GeV)] = 0.0038 ± 0.0006 GeV , U = 0.093 GeV . 



771^ = 0. 139 GeV, 



m c (2 GeV) 
m q (2 GeV) 



= 305 ± 59 , r b = 



m b (2 GeV) 
m q (2 GeV) 



1229 ±210, 



(36) 

(37) 



where the heavy quark mass at 2 GeV is determined using (32). The mixed condensate is similarly defined as 

m(qcrGq) = Mq m(qq) , 
where Mq = (0.8 ± 0.1) GeV [65] and m(qq) is as defined in (35 1. The gluon condensate is taken to be 

(aG 2 ) = (7.5 ± 2.0) x 10~ 2 GeV 4 [66] . 
Finally, the dimension-six quark condensate is 

a(qq) 2 = (5.8 ± 0.9) x 10~ 4 GeV 6 , (40) 

which implicitly includes deviation from ideal vacuum saturation [67 . In condensate contributions there are addi- 
tional factors of the quark mass that are not included in the definitions ( 35 ) or ( 38 ) , such as the factors of m 2k 



(38) 



(39) 



in (17 1, for instance. We define these masses in terms of the pole mass following the approach of Ref. J68], utilizing 



the known relationship between the pole mass and MS mass [551 E3 E3 IZ2] : 



to = m (p) 



1 



-log 



r 9 " 

m 

L 2 _ 


\ a(p)~ 



(41) 



where to (/j) and a (p) are determined via (32) and to is the one-loop MS charm or bottom quark mass. 

In order to extract mass prediction for heavy-light diquarks using (31 ) we must first establish a permissible range 



of values for the Borel scale r and the continuum threshold sq. We adopt the approach developed in Ref. 73 , 
whereby the Holder inequalities [7H [75] 



■J i. 



f(t)g(t)d[i 



< 



t 2 -, i/p 

\f(t)\ p dn 



Utt 



t 2 -, 1/9 

\g(t)\ q dn 



1 
P 



= 1, p,q> !• 



(42) 



are used to constrain the values of r and sq- The key observation of Ref. [73] is that because Imll (Q 2 ) is related to a 
physical hadronic spectral function via duality, Imll (Q 2 ) must be positive and hence it can serve as the integration 
measure in ( 42 ) . It can be shown that the sum rules ([28]) must satisfy 



7^2 (t, s ) /Hi (t,s ) 
fci (t , s ) /Tl (t , s ) 



> 1. 



gj (t ,_sq) /U 2 {t ,s q ) 
T^i {t > so) /Tli {r , so) 



> 1 



(43) 



where the first and second inequalities come from requiring that H (r ,sq) and Hi (r ,Sq) satisfy the Holder inequal- 
ities, respectively. The inequalities in (43) can be used to set a lower bound on the Borel mass Mb = 1/%/t" or to set 
a lower bound on the continuum threshold sq. The constraints set by the first inequality in (43) are more restrictive 



than those set by the second, hence we rely solely upon the first. We fix an upper bound on Mb by requiring that 
continuum contributions are less than 50% of total contributions to the sum rule 161 



/cont (T, Sq) 



fci (t, s ) /Tip (t, 8q) 
Hi (r, oo ) /H (r, oo ) 



(44) 



and require that /cont > 0.5. Using (43) and (44) we can define a range of Mb values over which the sum rule is 



considered reliable, i.e., the sum rule window. We also require that the mass prediction M (t , so) extracted from (31 ) 
exhibits r stability, that is, 



d_ 



M(r,s o )=0 



(45) 
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within the sum rule window. However, note that the bounds on the Borel scale that are determined using (43) 
and (44) will vary depending on the value of sq. Typically the sum rule window widens as sq is increased, 
we first seek a minimum value s, 



Thus 



„min 

S ! 



which we take to be the smallest value of sq in whose sum rule window r 
stability ( 45 ) is satisfied. If there are no values of Sq that exhibit r stability, we consider the sum rule to be unstable. 
Once the minimum value of sq has been found, we determine the optimal value Sg P using the following criterion: 



X 2 (so) 




so > 4" 



(46) 



The optimal value s is that for which (46 ) is minimized. We adopt a conservative approach where (46 1 is calculated 



over the sum rule window corresponding to the minimal value of sq. In some cases we can obtain an upper bound 



on the mass prediction (31 1 by taking so — > oo, however in order to extract such a bound the requirements described 
above must be satisfied. 



[Qq] J p M(GcV) M max (GeV) sf n (GeV 2 ) 



Mf in (GeV) 



A/£ iax 



(GeV) 



°Pt 



GeV 2 ) 



N 


0+ 


1.86 ±0.05 


2.02 


5.0 


1.2 


1.6 


5.0 




i+ 


1.87 ±0.10 


2.07 


5.0 


1.3 


1.6 


5.0 


N 


0+ 


5.08 ±0.04 


5.32 


30 


2.1 


3.8 


30 




1+ 


5.08 ±0.04 


5.32 


30 


2.2 


3.8 


30 



Table 4: Mass predictions and sum rule parameters for charm- light ([cq]) and bottom-light ([bq]) 
diquarks with positive parity. M max is an upper bound on the mass, determined from so — > oo. 
The minimal value of the continuum threshold is s™ 111 , the optimal value determined by (46) is 



„°P t 

b > 



We also determine the uncertainty in our mass predictions due to uncertainties in the QCD parameters. In order 
of significance these are r c 

and a(qq) 2 are insignificant in 



(37), m c (33), a(M T ) p3| , and M 2 (38) in the charm analysis, whereas in the bottom 
analysis they are r& (37), mj, (33), M<f(38), and a (Mz) (33). Uncertainties in (aG 2 ) 



both cases and we have made no attempt to estimate contributions to these uncertainties from higher loop effects. 
The resulting mass predictions and uncertainties for heavy-light diquarks with positive parity are summarized in 
Table [4] None of the negative parity heavy-light diquark sum rules exhibit r stability, therefore we have been unable 
to extract mass predictions in these channels. Fig. [6] shows the mass predictions for J p = + and 1 + charm-light 
diquarks while those for bottom-light diquarks are shown in Fig. [7J 



'«,/«o (GeV) 


2.15 




2.10 




2.05 


\ 


2.00 




1.95 




1.90 




1.85 







9 1.0 



1.1 



1.2 1.3 1.4 1.5 



^,/Ro (GeV) 



1.6 1.7 

tP — n+ 



M B (GeV) 




0.9 1.0 



1.1 1.2 1.3 1.4 



M B (GeV) 



Figure 6: Mass predictions for J = (left) and 1 (right) charm-light diquarks. Solid lines 
correspond to Sq P * , yielding the results in Table |4J In both plots the uppermost dashed lines 
correspond to so — > oo which provides an upper mass bound and the middle dotted line corresponds 
to s = 6.0 GeV 2 . 



11 



yjKj /Ko (GeV) 



-^K./Ro (GeV) 




M B (GeV) 




M„(GeV) 



Figure 7: Mass predictions for J p — + (left) and 1 + (right) bottom-light diquarks. Solid lines 



correspond to s 
correspond to so 



yielding the results in Table |4J In both plots the uppermost dashed lines 
oo and the middle dotted line corresponds to so = 35 GeV 2 . 



4 Conclusions 

In this paper we have used QCD Laplace sum rules to study heavy-light diquarks with J p — 0,1. Our cal- 
culations extend previous sum rule work |33j by including higher-loop perturbative contributions which necessitate 
renormalization of the diquark currents. We have successfully extracted mass predictions for positive parity charm- 
light and bottom-light diquarks, which are summarized in Table [4] However, the sum rules derived for negative 
parity channels are poorly behaved, and do not permit unambiguous mass predictions, similar to what was found for 
light diquarks [3"T] . 

The mass predictions for the J p = + and 1 + heavy-light diquarks are degenerate within uncertainty, as would 
be expected by heavy-quark symmetry |17j . Our predicted J p = + and 1 + charm- light diquark masses of 1.86 ± 
0.05 GeV and 1.87 ±0.10 GeV are in superb agreement with the constituent charm-light diquark mass of 1.93 GeV 
determined by Maiani et al. |T7j from a fit to the A (3872). Additionally, we predict both the J p = + and 
1 + bottom-light diquark masses to be 5.08 ± 0.04 GeV in reasonable agreement with the constituent bottom-light 
diquark mass of 5.20 GeV determined by Ali et al. [34] from a fit to the Yb(10890). Given the agreement between these 
constituent diquark masses and our QCD-based calculations, our results provide QCD support for the identification 
of the X(3872) and Fb(10890) as J PC = 1 ++ tctraquarks composed of diquark clusters. Furthermore, because 
the constituent heavy-light diquark is such an important input for constituent diquark models of tetraquarks, we 
interpret this agreement as indirect support for the predictions of these models. Specifically, our results strengthen 
the case for the tetraquark interpretation of the charged XYZ states -2^(3895), Z^(10610) and ^(10650). 

In this work we have focused on heavy-light diquarks so as to examine the constituent diquark masses determined 
in Refs. (TTJ |M]. However, the methods used in this paper could be extended to doubly- heavy diquarks to study 
diquark clustering within other tetraquarks or within heavy baryons. 
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